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A eultlvarlable control design procedure based on 
output feedback regulator (OFR) theory Is applied 
to the FIDO turbofan engine. Results for the OFR 
design are coasted to a design based on linear 
quadratic regulator (LQR) theory. This LQR design 
was (Attained as part of the FlOO Multivariable 
2 Cffittrol Synthesis (MFCS) program. In the MVCS 
jK progm the LQR feedback control was designed in a 
reduced diasnslon state space and then applied to 
the original system. However, the OFR feedback 
ctmtrol is designed In the full order state space 
and thus eliminates any need for model reduction 
techniques. Using the performance sieasure and 
control s^’ructure of the MVCS program LQR design, 
m eqvivalent OFR feedback control Is obtained. 

The flexibility of the OFR as a control design 
procedure Is demonstrated and differing feedback 
control structures are evaluated. 

INTRDIHJCTION 

Increased performance requirements, more com- 
plex engine configurations, and the possibility of 
onboard digital engine controllers have generated 
interest In advanced, mulclvarlable, engine con- 
trol system. To design these advanced control 
system engine manufacturers, as well as NASA and 
the Air Force, are Investigating multivariable 
design techniques (refs. 1 to 5). One program, 
the FlOO Multivariable Control Synthesis (MVCS) 
progrm, was Initiated to study the applicability 
of linear, quadratic regulator (LQR) theory to the 
design of FlOO turbofan engine controls (ref. 6). 
One aspect of the program Included obtaining con- 
stant feedback gains for regulation of engine 
ateady-state conditions. 

Engine data In the form of a nonlinear digi- 
tal simulation were used to generate linear models, 
each valid for small perturbations about an opera- 
tion point. However, the states corresponding to 
the linear models are not all physically measur- 
able for use in the feedback control law as would 
be required by LQR theory. Estimation of the un- 
masurable states from a masurable subset was 
judged, within the program hardware constraints, 
to be too coi^lex to Inplemnt. Thus, the control 
design procedure selected first re&iced via doml- 
V nant suds approximation techniques (ref. 6) the 
order of the mdel state space to Include only a 
masurable stbspace. Then a state feedback regu- 
lator was designed using LQR theory In the reduced 
state space. This procedure gave an implementable 


control design. Moreover, a real-tlm hybrid com- 
puter simulation evaluation (ref. 7; has demon- 
strated (1) the ability of the resultant control 
to perform the prescribed control function and (2) 
the flexibility of the LQR design process for de- 
termining engine controls. 

However, the model reduction process, or more 
specifically the selection of the masurable sub- 
space, was not straightforward and complicated the 
overall design procedure. In particular the selec- 
tion of the reduced state space defines a trade-off 
between control performance and the control liq>le- 
mntatlon complexity. A definitive answer to the 
adequacy of this trade-off would require many Iter- 
ations through the design procedure and, therefore, 
Biany model reductions. Additionally, there Is no 
a priori assurance that any approximation technique 
selected will yield good (or even acceptable) re- . 
suits when the resultant control design Is applied 
to the original system. Therefore, It la desirable 
to determine a design process that would incorpor- 
ate the utility of Che LQR theory but would elimi- 
nate the model reduction step. 

The purpose of this paper Is to propose such 
a design process based on output feedback regu- 
lator (OFR) theory (ref 9) and to Investigate the 
acceptability of this process for designing engine 
controls. The OFR design process incorporates a 
linear, quadratic formulatlcn but would use the 
full order linear model directly to determine con- 
stant feedback gains for an appropriate output vec- 
tor (som vector In the measurable subspace). 

Thus, Che model reduction step Is eliminated. 
Additionally, Che OFR formulation allows an effi- 
cient study of the performance trade-off liqportant 
In any realistic design. These benefits are not 
obtained without seme cost, however. Solutions to 
the OFR problem satisfy only a necessary optimality 
condition and therefore are suboptlmal «dien com- 
pared to the necessary and sufficient conditions of 
the full state feedback LQR forimilatlon. Recall, 
however, that the LQR process as applied in the 
MVCS program used a reduced order state space and, 
therefore. Is also suboptlmal with respect to the 
original system. Thus, OFR suboptlmallcy la not 
regarded as too significant a problem. More Impor- 
tantly, OFR solutions are not guaranteed to exist 
for any arbitrary output vector. Thus, OFR solu- 
tions must be shown to exist for the desired con- 
trol structure before proceeding with the design. 

The OFR design process was Investigated by 
applying It to a 17th order linear, operating point 
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■odal of Cha FlOO turbofan engine. Tala model 
repreaenta an engine Intermediate power condition 
In a aea level, static environment. First, a 
brlaf ravlew of OFR theory la given. Then, com- 
parisons are made between the resultant OFR feed- 
back control and the MVCS program LQR control de- 
signed at the same condition. Finally, the flexi- 
bility of the OFR approach Is demonstrated by 
evaluating the Improvement In performance when 
additional engine variables are measured and In- 
cluded In the control configuration. 

HATUEMATICAL DEVELOPMENT 
Given a time Invariant linear system 


X ■ Ax + Bu, x(0) “ Xq 


( 1 ) 


y •. Cx 

where the Initial 'jtate, Xq, Is a zero-mean random 
variable with covariance, Xq. The OFR problem Is 
to find the time Invariant feedback law 


u - -Fy (2) 

which minimizes 


environment. The lapi:t and full state vectors 
that represent the FlOO engine variables are given 
In Table 1. The numerical description of this 
model la given In reference 6. For comparison pur- 
poses the same performance weightings used for the 
MVCS program LQR design at this condition were used 
In the OFR design. These numerical weightings are 
also given in reference 6. 

The LQR and OFR procedures used the same In- 
put vector but the LQR approach selected a reduced 
state vector as shown by equation (6). 



if we let the output vector y be the reduced 
order state vector 


y ■ X “ Cx 


(9) 


with 


J 



(x^Qx + u^Ru)dt 


Q 2. 0 
R > 0 


(3) 


(4) 


Necessary conditions for optimality and computa- 
tional solutions are derived for this problem 
statement In references 8 and 9 and for Its dis- 
crete time counterpart In reference 4. The neces- 
sary conditions are 


then the regulator structure when applied to the 
original system Is the same for both the OFR and 
LQR designs. Computationally, the feedback matrix, 
F, of equation (2), was solved by approximating 
the continuous system of (1) by its discrete 
counterpart and applying the numerical algorithms 
of reference 4. The sampling period was selected 
as T •• 0.0001 secs/cycle to give a close approx- 
imation to the continuous system since the open 
loop eigenvalues are greater than -600 rad/sec. 

With this approximation the output feedback matrix 
could be directly compared to that obtained via 
the LQR design. Al?o, the Initial condition covar- 
iance matrix was somewhat arbitrarily set as 
Xq • I for all the OFR designs of this paper. 


KAjj + aJk + Q + C^F^RFC - 0 

(5a) 

LaJ + .ij)L + Xq - 0 

(5b) 

F - r"^B^KLC^(CLC^)“^ 

(5c) 

where 

Ajj - A - RFC 

(6) 

The solution of (5) gives F, and the 
value of J for this F is given by 

suboptlmal 


Tne numerical procedure is Iterative in nature 
and requl ■•.s a stabilizing output feedback matrix 
as an Initial guess. Sirce the FlOO linear models 
are open loop stable, an Initial guess of zero for 
each feedback element can be conveniently selected 
to start the Iterative process. With this Initial 
guess the output feedback matrix, F, was found for 
this model quite easily. The feedback matrices 
and the closed loop eigenvalues for the OFR and 
LQR approaches are compared In Table II. Also, 
Included are the final cost function values. The 
results compare quite favorably and It can be seen 
that control performance equivalent to the LQR de- 
sign could have been directly determined by the 
OFR design process. 


J - j Tr(KXg) 


(7) 


OFR DESIGN FLEXIBILITY 


OFR DESIGN FOR FlOO 

r^e OFR design procedure outlined in the pre- 
ceding section was applied to the 17th order 
linear time Invariant model of engine operation 
point, small perturbation, dynamics for an Inter- 
mediate power condition In a sea level, static 


T!ie OFR design procedure offers the usual 
features of LQR design, plus the additional capa- 
bility of Investigating different structures (via 
the C matrix) without altering the basic dynamic 
model of the system (A and B) . To demonstrate this 
flexibility the FlOO model was studied to determine 
If different combinations ot measurable engine 


variable* could result In better control perform- 
ance or a simplified control structure, given the 
previously defined performance weightings. With 
the performance measure of (7) It Is possible to 
compare control effectiveness for different feed- 
back control structures. 

Three different structures were compared for 
the given engine model. The output vectors for 
these structures are 

Structure #1 

y - X (10) 

Structure 12 



SUMMARY 

This paper has presented a multivariable con- 
trol design procedure based on OFR theory that can 
be used to design operating point controls for Jet 
engines. The procedure utilizes the benefits of a 
linear, quadratic approach but eliminates the 
model reduction step required by LQR theory when 
the number of measurable engine variables Is less 
than Che dimension of the linear model state space, 
.tddlclonally , sensor failure and complexity- 
performance trade-off studies. In the form of dif- 
fering feedback structures, can be handled quite 
readily with this formulation. Results were ob- 
tained by applying the OFR procedure to a linear 
17th order model of the FIDO Curbofan engine and 
compared Co Chose obtained by a reduced dimension 
LQR approach. 
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TABLE 1. DEFINITIONS OF CONTROL AND STATE VECTORS 
FOR FIDO MNEAR HODEL 


“1 


P * 

CoBUBsnded fuel flow 

“2 


Nozzle Jet 

“3 

- 

Inlet guide vane position 

“4 


High compressor variable stator position 

“5 

• m 


Compressor bleed flow 

^ • 

■*r 


'Fan speed 

*2 


Compressor speed 

*3 


Compressor discharge pressure 

*4 


Inter-turbine volume pressure 

*5 


Augmentor pressure 

*6 


Fan inner diameter discharge temperature 

*7 


Duct temperature 

*8 


Compressor discharge temperature 

*9 


Fast response burner exit temperature 

*10 


Slow response burner exit temperature 

*11 


Burner exit total temperature 

*12 


Fast response fan turbine inlet temperature 

*13 


Slow response fan turbine inlet temperature 

*14 


Fan turbine exit temperature 

*15 


Duct exit temperature 

*16 


Augmentor exit temperature 

.*17 


Main bur.ner fuel fl.'«w 


1 



TABLE II. COMPARISON OF LQR AND OFR DESIGNS 
OFR Feedback Matrix 


“1 

“2 

U3 

“4 

“5 


0.336560 

.580896 

35.5966 

1.53157 

-0.102381 


1.42533 

.191406 

3.02362 

-66.1373 

-0.803223 


0.284484 

-0.332299 

-3.22347 

12.7661 

-1.98148 


-0.244595-02 

-0.528128 

-8.75157 

1.37397 

-0.127121 


‘17 

0.196502 

-0.526202-01 

1.28023 

-2.27481 

-0.735746-01 


LQR Feedback Matrix 



“1 

“2 

“3 

“4 

“5 


0.497890 

.690374 

36.1351 

1.18785 

.837923-01 

OFR Coat (J) 


2.42002 

.211490 

-6.78935 

-23.6908 

-0.992363 


-0.101360 

.421472-01 

9.98711 

21.8884 

-0.946234 


0.120072 

-0.553484 

-8.68853 

2.04405 

-0.258270 


‘17 

0.240518 

-0.765097-01 

.547046 

.216411 

-0.120667 


LQR Coat (J) 


The optlaaia coet 
55127.5 


Trace (Rlcattl matrix) 



OFR Cloeed Loop Elgenvaluea 


*5 


*17 


-554.703 

+ J* 

15.3219 

*1 

-349.868 

+ 

-554.703 

+ j! 

-15.3219 

*2 

-10.0215 

+ 

-52.7604 

+ J 

3.14720 

*5 

-10.0215 

+ 

-52.7604 

+ j! 

-3.14720 

*17 

-10.6641 

+ 

-46.5396 

+ J 

.000000 

*1 

-6.13386 

+ 

-40.2390 

+ J* 

.000000 

1 



-24.6247 

+ J.. 

.000000 




-18.1467 

+ J* 

3.89611 




-18.1467 

+ J 

-3.89611 




-21.3023 

+ J* 

.000000 




-18.8500 Q 

+ J* 

.000000 




-9.65625 

+ j! 

2.31956 




-9.65625 

+ JV 

-2.31956 




-0.669755 

+ J 

.000000 




-6.14437 

+ J* 

.000000 




-1.97937 

+ J* 

.000000 




-10.5293 

+ J* 

.000000 





The optimum coat " Trace (Rlcattl matrix) 
55128.8 

LQR Cloaed Loop Elgenvaluea 


0.000000 

2.96296 

-2.96296 

.000000 

.000000 




TABLE 111. COMPARISON OF THREE CONTROL STRUCTURES 




Structure II 


*1 

*2 

*3 

*5 

*17 

U.33a03 

1.42786 

0.241943 

0.106119-02 

0.196665 

.5fi3311 

. 183079 

-0.278551 

-0.531235 

-0.538824-01 

35.5320 

2.94085 

-2.64476 

-8.81418 

1.27011 

1.28780 

-65.6406 

11.3587 

1.54749 

-2.22466 

-0.800610-01 

-0.846198 

-1.73410 

-0.153524 

-0.784924-01 


Structure 12 


*1 


*3 


*5 


X 


17 


“1 

“2 

U3 

“5 


0.724593 3.35082 

.621300 .672284-01 

35.7266 1.91687 

-13.7957 -124.472 

-0.196554 -2.82217 


-0.574654 -0.129167 
-0.603105 -0.973783-01 
-9.97863 .477242 

27.1895 13.8395 

.122541 .186423 


Structure 13 


*1 

*2 

*3 

*5 

*17 

*14 

0.324920 

1.41944 

0.251627 

-0.676767-02 

0.184383 

0.548166-01 

.582912 

.182919 

-0.278377 

-0.531972 

-0.546010-01 

.617017-02 

35.4674 

2.84734 

-2.58472 

-8.93478 

1.11033 

.806580 

1.50482 

-65.4595 

11.2330 

1.80488 

-1.90962 

-1.65499 

-0.795668-01 

-0.845408 

-1.73433 

-0.152529 

-0.770434-01 

-0.655999-02 


OFR Cost (J) 


Structure 1 55 237.8 
Structure 2 58 164.8 
Structure 55 233.8 
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